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ABSTRACT 


We present the observational results of the y-ray blazar, 3C 66A, at 2.3, 8.4, and 
22 GHz at 4 epochs during 2004-05 with the VLBA. The resulting images show an 
overall core-jet structure extending roughly to the south with two intermediate breaks 
occurring in the region near the core. By model-fitting to the visibility data, the north- 
most component, which is also the brightest, is identified as the core according to its 
relatively flat spectrum and its compactness. As combined with some previous results 
to investigate the proper motions of the jet components, it is found the ki nematics of 
3C 66A is quite complicated with components of inward and outward, subluminal and 
superluminal motions all detected in the radio structure. The superluminal motions 
indicate strong Doppler boosting exists in the jet. The apparent inward motions of the 
innermost components last for at least 10 years and could not be caused by new-born 
components. The possible reason could be non-stationarity of the core due to opacity 
change. 


Subject headings: radio continuum: galaxies— galaxies: individual (3C 66A)— galaxies: jet 


-3- 


1. Introduction 


The radio source, 3C 66A (also known as B0219+428, NRAO102, 4C 42.07), is a 
well-known blazar, which is further sub-classified as an intermediate synchrotron peaked BL Lac 
object (IBL) ( Abdo et al. 2010). Since its optical counterpart was identified by Wills & Wills 
(1974), 3C 66A was studied over a wide frequency range, from radio to y-ray. It is also one of 
the sources detected in TeV (Aliu et al. 2009, Acciari et al. 2009). Like other blazars, 3C 66A 
exhibits violent variabilities at radio, IR, optical, and X-ray (e.g. Bottcher et al. 2005 hereafter 
B05), and in recent years, flares at y -ray wavelengths have also been detected (e.g. Swordy 2008, 
Abdo et al. 2011). 


On kpc-scales, 3C 66A shows a 6 arcsecond extended structure at the position angle (P.A.) ~ 
170° (Price et al. 1993) with two weak lobes, one at 15" south of the central region, the other at 
10” with P.A. about -20° ( Taylor et al. 1996). Very Long Baseline Interferometry (VLBI) images 
from 2.3 to 43 GHz show a typical core-jet structure on parsec scales (e.g. Taylor er al. 1996, 
Jorstad et al. 2001, Marscher et al. 2002, Cai et al. 2007, hereafter C07). The kinematics of 3C 
66A jet has been studied via multi-epoch VLBI observations with superluminal motions detected 
for several jet components (e.g. Jorstad et al. 2005, Britzen et al. 2008). 


The redshift value of 3C 66A is adopted to be z=0.444(Miller et al. (119781) ) by several 


previous works (e.g. B05 and C07) . This value was measured based on the detection of only one 
single line and is argued to be quite uncertain (e.g. Bramel et al. 2005). Recently, several redshift 
ranges for this source were obtained by different methods, e.g. z < 0.58 by Yang & Wang (2010), 
z > 0.30 by Shaw et al. (2013), 0.335 < z < 0.41 at 99 percent confidence and 0.335 < z < 0.444 
at 99.9 percent confidence by Fumiss et al. (2013), and z > 0.42 by Stadnik & Romani (2014). In 
this work, we adopt 0.335 < z < 0.444 as the redshift range of 3C 66A. 


In this paper we present our multi-frequency, multi-epoch VLBA observation results of 
3C 66A. The observations and data reduction are described in Sect. 2. The results, including 
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the images, spectra, and proper motion analysis are presented in Sect. 3. In Sect. 4, a detailed 
discussion is made with a summary followed in Sect. 5. In this paper,we use the ACDM 
cosmological model with H 0 = 71 km s _1 Mpc -1 , = 0.27, and Q A = 0.73 (Spergel et al. 
2003). 


2. Observations and data reduction 

3C 66A was observed as phase-reference calibrator for the VLBA projects BS087 and 
BS144, which aimed at searching for evidence of a super-massive black hole binary in the nearby 
radio galaxy 3C 66B (Sudou et al. 2003, Zhao et al. 2011). These projects were carried out during 
2001-2002 (Session 1) and 2004-2005 (Session 2), respectively. Session 1 includes 6 observations 
at 2.3 and 8.4 GHz, and 2 observations at 22 GHz at the first 2 epochs. The results of Session 1 
were published in C07. During Session 2, 3C 66A was successfully observed at 2.3, 8.4, and 22 
GHz at all 4 epochs, 2004.80, 2005.05, 2005.35, and 2005.54. The total on-source time of 3C 
66A at each epoch is about 80 min at 2.3 GHz, 30 min at 8.4 GHz, and 40 min at 22 GHz. The 
data were recorded in VLBA format with 2 intermediate frequency (IF) bands at 2.3 and 8.4 GHz 
and 4 IF bands at 22 GHz. Each IF has a bandwidth of 8 MHz. The total recording rate was 128 
Mbps using a 2-bit sampling mode. The observed data were correlated with the VLBA correlator 
in Socorro. 

The output data were then loaded into the NRAO AIPS package for further processing. 
Amplitude calibration was performed by using the measured system temperature and antenna 
gain curves. Phase calibrations were made in a standard way given in AIPS COOKBOOK, except 
that additional correction for all of the observations in Session 2 was made to correct the possibly 
inaccurate Earth orientation parameters (EOPs) used by the VLBA correlator from 2004 to 2005. 
The residual phase delays and delay rates that were mainly due to antenna position errors and 
atmospheric phase fluctuations were calibrated using global fringe-fitting (Schwab & Cotton 
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1983). Bandpass calibration was performed for higher sensitivity using 0133+476 for all epochs. 
The output data from AIPS were then read into the Caltech DIFMAP package for hybrid mapping 
and model fitting. Circular Gaussian models were used in model fitting, which may reduce the 
number of free parameters and give more consistent centroid positions across the epochs (Lister 
et al. 2009). The corresponding er rors fo r the fitted parameters were estimated by using the 


expressions described in Lee et al. (120081) . except the flux density errors were taken to be 10% at 


2.3 and 8.4 GHz, and 20% at 22 GHz, and the errors for the core separations of the components 
were taken to be at least 0.8 mas at 2.3 GHz, 0.2 mas at 8.4 GHz, and 0.06 mas at 22 GHz. Such 
criterions were adopted based on some previous works on other radio sources (e.g. C07, Lu et al. 
2011 ). 


3. Results 

In Fig. [H we show the resulting contour images of 3C 66A at all 4 epochs at 2.3, 8.4, and 
22 GHz, with specific parameters shown in Table. [TJ An overall core-jet emission distribution 
can easily be found at low frequency of 2.3 GHz with higher dynamic range and larger emission 
structure, while more specific structures are obtained at 22 GHz with higher resolution. The radio 
structure at all these frequencies and epoches shows a consistent intensity distribution with a 
compact emission feature lying at the extreme end of a onesided jet, and others extending to the 
south. From the 8.4 GHz radio structures at all 4 epochs, a prominent bending was found at ~ 4 
mas. And in the 22 GHz images, a minor bending occurs at ~ 1.2 mas. 

To better understand its intensity distribution down the jet, model fitting was done to the 
remarkable jet features by using the modelfit task in DIFMAP. All the resulting components are 
marked with circular signs with their sizes scaling to their true emission regions. Table. [2] shows 
the designations of each jet component, the corresponding flux densities and uncertainties at all 
epochs and frequencies. The designation of each component mainly follows the convention in C07 



- 6 - 


Table 1: Description of VLB A maps of 3C 66 A shown in Fig. Q] 





Restoring Beam 


V 

Epoch 

Speak 

Major 

Minor 

PA. 

Contours 

(GHz) 

(yr) 

(Jy/beam) 

(mas) 

(mas) 

(deg) 

(mJy/beam) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

2.3 GHz 

2004.80 

0.453 

5.86 

3.77 

-0.74 

0.9x(-1,1,2,4,8,16,32,64,128,256,512) 


2005.05 

0.476 

5.89 

3.58 

-3.21 

0.9x(-l,l,2,4,8,16,32,64,128,256,512) 


2005.35 

0.424 

7.51 

4.25 

-18.5 

1.0x(-1,1,2,4,8,16,32,64,128,256) 


2005.54 

0.465 

6.31 

3.96 

-0.72 

0.9x(-l,l,2,4,8,16,32,64,128,256,512) 

8.4 GHz 

2004.80 

0.451 

1.61 

1.10 

5.22 

1.2x(-1,1,2,4,8,16,32,64,128,256) 


2005.05 

0.466 

1.51 

1.01 

0.472 

1.0x(-1,1,2,4,8,16,32,64,128,256) 


2005.35 

0.496 

2.11 

1.18 

-18.2 

1.4x(-1,1,2,4,8,16,32,64,128,256) 


2005.54 

0.433 

1.56 

1.04 

-1.05 

1.2x(-1,1,2,4,8,16,32,64,128,256) 

22.2 GHz 

2004.80 

0.464 

0.829 

0.575 

-0.43 

2.0x(-1,1,2,4,8,16,32,64,128) 


2005.05 

0.421 

0.614 

0.371 

-3.29 

1.5x(-1,1,2,4,8,16,32,64,128,256) 


2005.35 

0.488 

0.991 

0.761 

3.72 

2. lx(-1,1,2,4,8,16,32,64,128) 


2005.54 

0.367 

0.629 

0.378 

-6.29 

1.8x(-1,1,2,4,8,16,32,64,128) 


Notes:(l) Observing frequency; (2) Observing epoch; (3) Peak flux density; (4), (5), (6) Parame¬ 
ters of the restoring Gaussian beam: the full width at half maximum (FWHM) of the major and 
minor axes and the position angle (PA.) of the major axis. (7) Contour levels of the map. The 
lowest contour level is three times the rms noise in the map. 
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for the convenience of proper motion analysis. The counterpart for each component at different 
frequencies are mainly determined by their relative positions to the core. A newly identified 
component at 8.4 and 22 GHz which lies at about 4 mas away from the northmost component is 
labeled as n. This component was not detected in the 2001-2002 session 1 observations of C07, 
which may be due to its weakness and relative closeness to component d during the observing 
campaign. 


3.1. Spectrum and identification of components 

The detection of components k, b, d, e, and/ at more than one frequency allows us to make 
spectral analysis of these components. Fig. [2] (top) displays the spectra of these components 
at the first epoch (2004.80) based on the flux densities obtained by model-fitting of Gaussian 
components to the visibility data. The spectral indices (5 oc y“) at all epochs are shown in Fig. [2] 
(bottom), among which the spectral indices of components k and d are obtained through fits to 
the data at all the 3 frequencies, and the others are obtained using data at 2 frequencies at which 
the components were detected. The component k at the north end of the whole radio structure 
shows the flattest spectrum {a ~ 0) among all these components, which enables us to identify this 
component as the core. By contrast, the other components show relatively steep spectra {a < 0). 

We also estimate the brightness temperatures of all the components in the source frame by 
using the expression (Shen et al. 1997) 

T h = 1.22 x 10 12 5 v" 2 d~ 2 (1 +z)K , (1) 

where S is the flux density in Jy, v is the observational frequency in GHz, d is the diameter in 
mas, and z is the redshift. The brightness temperatures were shown in the last collum in Table. 13 
Component k has the highest brightness temperature among all components, which greatly 
exceeds the Compton limit of 10 12 AT and hence gives further supports to our identification of the 


core component, as is shown in Table El 


3.2. Proper motions of the jet components 

C07 studied the ki nematics of the jet components in 3C 66A at 2.3 and 8.4 GHz using 
Session 1 data and no obvious proper motion was detected within the uncertainties because of the 
short time coverage (~1.2yr). They combined their results with some previous works (Jorstad 
et al. 2001, Jorstad et al. 2005, and B05) and they found the kinematics of 3C 66A is rather 
complicated, yet the observations were carried out at different frequencies. By combining the data 
of the two sessions, we are able to study the ki nematics of 3C 66A over a much longer time range 
(more than 4 years). The 22 GHz preliminary results were already published in the proceeding 
of IAU Symposium 290 (Zhao et al. 2013). In Fig. |3j we show the plots of separations of jet 
components from the core that have measurements at both sessions as a function of time. The 
apparent velocity obtained by linear fit to the separations are shown in Table. |3] We compared our 
22 GHz results with that of B05 which was obtained in the time gap between our two Sessions 
(2003-2004). The model fitted positions of components a, b, and d are in good agreement with 
C3, C2, and B2, respectively, in their work (Fig. [3]c), while component c seems to be divided into 
two components (B3 or Cl). So we also included B05 results (Comp. C3, C2, & B2) in our 22 
GHz proper motion study. 

Our results show that the speeds of components in the jet of 3C 66A are much slower than 
previously reported (e.g. Jorstad et al. 2001, 2005). The results confirmed that the kinematics of 
this source is complicated, i.e. some components show superluminal motions, some show inward 
motions, and the others are nearly stable. As shown in Fig. [3 there are no obvious proper motions 
found for component/ at > lcr level. Components a, b, and g show inward motions at > lcr level 
while components c, d, and e move downstream the jet superluminally at > 3cr level. 
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Table 2. Parameters of the Gaussian model fit components of the images shown in Fig. [Q 


Freq. (GHz) Comp. 

S(Jy) 

R (mas) 

6 (deg) 

Major (mas) 

Tb (K) 

Epoch 2004.80 

S 

k 

0.362+0.036 

0 

0 

0.35+0.00 

1.01 x 10 12 


d 

0.127+0.013 

2.26+0.80 

-176.9+0.1 

0.82+0.01 

6.28 x 10 10 


e 

0.049+0.005 

7.08+0.80 

170.1+0.5 

2.94+0.13 

1.90 x 10 09 


f 

0.061+0.006 

13.03+0.80 

162.9+1.0 

4.92+0.47 

8.45 x 10° 8 


g 

0.032+0.003 

21.49+0.86 

167.4+2.3 

9.24+1.72 

1.23 x 10° 8 

X 

k 

0.401+0.040 

0 

0 

0.09+0.00 

1.11 x 10 12 


b 

0.076+0.008 

0.68+0.20 

-171.7+0.6 

0.30+0.01 

2.19 x 10 10 


d 

0.088+0.009 

2.26+0.20 

-173.4+0.4 

0.58+0.03 

6.66 x 10° 9 


n 

0.018+0.002 

4.20+0.23 

179.8+3.2 

1.99+0.47 

1.15 x 10° 8 


e 

0.033+0.003 

8.01+0.28 

170.0+2.0 

2.82+0.56 

1.02 x 10° 8 


f 

0.037+0.004 

13.43+0.62 

164.4+2.7 

4.26+1.25 

5.06 x 10° 7 

K 

k 

0.402+0.080 

0 

0 

0.04+0.00 

9.02 x 10“ 


a 

0.077+0.015 

0.27+0.06 

-171.4+0.5 

0.13+0.00 

1.61 x 10 10 


b 

0.030+0.006 

0.85+0.06 

-174.7+0.8 

0.35+0.02 

8.63 x 10° 8 


c 

0.007+0.001 

1.72+0.06 

-166.0+1.3 

0.33+0.08 

2.21 x 10° 8 


d 

0.058+0.012 

2.35+0.06 

-172.8+0.4 

0.42+0.03 

1.19 x 10° 9 


n 

0.006+0.001 

3.99+0.15 

-174.0+2.2 

0.67+0.31 

4.91 x 10° 7 

Epoch 2005.05 

S 

k 

0.383+0.038 

0 

0 

0.13+0.00 

7.63 x 10 12 


d 

0.130+0.013 

2.32+0.80 

-175.7+0.1 

0.87+0.01 

5.65 x 10 10 
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Table 2—Continued 


Freq. (GHz) 

Comp. 

S(Jy) 

R (mas) 

6 (deg) 

Major (mas) 

Tb (K) 


e 

0.050+0.005 

7.06+0.80 

170.6+0.7 

2.97+0.18 

1.87 x 10 09 


f 

0.061+0.006 

12.82+0.80 

163.3+1.1 

4.89+0.49 

8.51 x 10° 8 


g 

0.034+0.003 

20.94+0.80 

166.9+2.2 

8.47+1.60 

1.57 x 10° 8 

X 

k 

0.430+0.043 

0 

0 

0.10+0.00 

1.05 x 10 12 


b 

0.063+0.006 

0.66+0.20 

-174.1+0.6 

0.32+0.01 

1.51 x 10 10 


d 

0.076+0.008 

2.30+0.20 

-173.9+0.4 

0.60+0.03 

5.17 x 10° 9 


n 

0.014+0.001 

3.96+0.20 

-178.5+2.5 

1.65+0.35 

1.27 x 10° 8 


e 

0.032+0.003 

7.58+0.32 

169.2+2.4 

2.87+0.64 

9.75 x 10° 7 


f 

0.039+0.004 

13.22+0.53 

163.4+2.3 

4.16+1.06 

5.64 x 10° 7 

K 

k 

0.378+0.076 

0 

0 

0.06+0.00 

3.48 x 10“ 


a 

0.078+0.016 

0.26+0.06 

-173.5+0.4 

0.14+0.00 

1.38 x 10 10 


b 

0.030+0.006 

0.73+0.06 

-174.8+0.6 

0.27+0.01 

1.55 x 10° 9 


c 

0.006+0.001 

1.50+0.06 

-170.5+0.6 

0.22+0.03 

4.41 x 10° 8 


d 

0.042+0.008 

2.42+0.06 

-173.0+0.7 

0.47+0.06 

7.06 x 10° 8 


n 

0.007+0.001 

3.67+0.31 

-173.9+4.8 

1.41+0.62 

1.29 x 10° 7 


Epoch 2005.35 


0.317+0.032 

0 

0 

0.05+0.00 

4.18 x 10 13 

0.126+0.013 

2.11+0.80 

-175.1+0.0 

0.76+0.00 

7.33 x 10 10 

0.051+0.005 

6.84+0.80 

170.9+0.5 

2.89+0.13 

2.02 x 10° 9 

0.062+0.006 

12.82+0.80 

163.1 + 1.1 

4.75+0.48 

9.21 x 10° 8 

0.034+0.003 

21.00+0.80 

166.5+2.2 

9.02+1.58 

1.38 x 10° 8 


g 
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Table 2—Continued 


Freq. (GHz) 

Comp. 

S(Jy) 

R (mas) 

6 (deg) 

Major (mas) 

Tb (K) 


X 

k 

0.427+0.043 

0 

0 

0.14+0.00 

5.79 x 10 11 


b 

0.093+0.009 

0.59+0.20 

-174.6+0.3 

0.39+0.01 

1.55 x 10 10 


d 

0.069+0.007 

2.29+0.20 

-174.2+0.2 

0.54+0.02 

6.01 x 10 09 


n 

0.013+0.001 

3.48+0.20 

-178.1+0.7 

0.91+0.08 

3.89 x 10° 8 


e 

0.040+0.004 

7.61+0.26 

169.8+2.0 

3.14+0.53 

1.00 x 10° 8 


f 

0.035+0.004 

13.25+0.39 

162.5 + 1.7 

3.50+0.77 

7.19 x 10° 7 

K 

k 

0.392+0.078 

0 

0 

0.02+0.00 

5.94 x 10 12 


a 

0.073+0.015 

0.21+0.06 

-163.2+0.2 

0.03+0.00 

2.59 x 10“ 


b 

0.066+0.013 

0.59+0.06 

-174.1+0.3 

0.15+0.01 

1.09 x 10 10 


c 

0.011+0.002 

1.67+0.06 

-171.5+0.1 

0.14+0.01 

2.16 x 10° 9 


d 

0.044+0.009 

2.45+0.06 

-174.5+0.4 

0.46+0.03 

7.67 x 10° 8 


n 

0.006+0.001 

3.69+0.13 

-179.4+2.1 

0.88+0.27 

2.84 x 10° 7 

Epoch 2005.54 

S 

k 

0.391+0.039 

0 

0 

0.30+0.00 

1.45 x 10 12 


d 

0.105+0.011 

2.38+0.80 

-175.5+0.2 

0.92+0.01 

4.14 x 10 10 


e 

0.045+0.005 

7.09+0.80 

170.9+0.6 

3.02+0.14 

1.65 x 10° 9 


f 

0.062+0.006 

12.71+0.80 

162.8+0.9 

4.64+0.41 

9.63 x 10° 8 


g 

0.035+0.003 

20.61+0.82 

166.9+2.3 

9.23+1.64 

1.35 x 10° 8 

X 

k 

0.388+0.039 

0 

0 

0.13+0.00 

5.40 x 10“ 


b 

0.075+0.008 

0.68+0.20 

-176.3+0.5 

0.24+0.01 

3.19 x 10 10 


d 

0.062+0.006 

2.30+0.20 

-174.7+0.5 

0.60+0.04 

4.35 x 10° 9 
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4. Discussion 
4.1. Morphology 

Our VLBI images of 3C 66A show similar morphology as previous works, which is a typical 
one-side jet with two bendings at about 1.2 and 4 mas from the core. Jet bendings or non-linear 
morphologies are very common in BL Lacs (e.g. Karouzos et al. 2012). Jorstad et al. (2005) 
suggest that the first bending in 3C 66A jet occurs smoothly from 1.2 mas up to 3 mas. The 
former works mentioned above all show that the flux densities of the components at 22 and 43 
GHz decrease along the jet but start to increase at 1.2 mas up to 2.6 mas. Such a coincidence 
suggests that the observed bending may be related to the cause of the local flux density increase 
in the jet, which points directly to the continues decrease of the viewing angles and hence the 
increase of Doppler boosting effect. The model fitting results of components k-d in our 22 GHz 
maps also support this argument. Bendings could also be originated by interaction with ambient 
medium (e.g. Alberdi et al. 2000), but this is less likely to be the case for the first bending in 3C 
66A because there are no stationary components in this area (see Sec. 13.21) . 

As measured in the other blazar objects, the detectable radio structure of 3C 66A gets smaller 
with frequency. The restored emission region at 22 GHz is almost 6 times smaller than that at 
2.3 GHz as a result of the much higher resolution plus the decreased surface brightness at 22 
GHz. Note that the brightness temperature drops with frequency, indicating its decreased surface 
brightness at higher frequency, as is suggested in Kellermann et al. (2004). 


4.2. Kinematics 

We fit the core separations of the jet components of 3C 66A for the two Sessions, which 
indicates the kinematics of 3C 66A is very complicated. We have detected apparent inward 
motions as well as superluminal motions. The superluminal moving components (c, d, e) are found 
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Table 2—Continued 


Freq. (GHz) 

Comp. 

S(Jy) 

R (mas) 

e (deg) 

Major (mas) 

Tb (K) 


n 

0.018+0.002 

3.70+0.20 

-177.9+2.7 

1.51+0.35 

1.94 x 10 08 


e 

0.032+0.003 

7.75+0.30 

170.5+2.2 

2.85+0.60 

9.96 x 10 07 


f 

0.037+0.004 

12.98+0.45 

161.8+2.0 

3.53+0.89 

7.44 x 10 07 

K 

k 

0.306+0.061 

0 

0 

0.04+0.00 

6.19 x 10“ 


a 

0.094+0.019 

0.24+0.06 

-164.9+0.4 

0.12+0.00 

2.31 x 10 10 


b 

0.043+0.009 

0.63+0.06 

-175.1+0.5 

0.21+0.01 

3.73 x 10 09 


c 

0.008+0.002 

1.58+0.06 

-172.5+2.0 

0.50+0.11 

1.21 x 10° 8 


d 

0.037+0.007 

2.50+0.06 

-173.9+0.5 

0.45+0.05 

6.62 x 10° 8 


n 

0.004+0.001 

3.36+0.06 

-179.5+0.3 

0.13+0.04 

8.45 x 10° 8 
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3C 66A at 2.3GHz 
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Fig. 1.— The naturally-weighted VLB A maps of 3C 66 A at 2.3 GHz (top), 8.4 GHz (middle), 
and 22.2 GHz (bottom). The circles superimposed on the maps represent the Gaussian model 
components listed in Table [2l whose names are labeled. The crosses represent the restoring beams 
in Table Q] 
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Fig. 2.— Top: Component spectra of 3C 66A on 2004.80. Bottom: spectral indices of jet compo¬ 
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Fig. 3.— Separations from the core as functions of time for the jet components in 3C 66A at 2.3 
GHz (top), 8.4 GHz (middle), and 22 GHz (bottom). The 2001-2002 results are taken from C07, 
2003-2004 results at 22 GHz are taken from B05. The open symbols represents the components 
that cannot be cross identified. 
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Table 3: The proper motions of the jet components in 3C 66A 



2.3 GHz 

8.4 GHz 

22.2 GHz 


masyr -1 

fiappiz = 0.444) 

V- 

masyr -1 

Pappiz = 0.444) 

mas yr _1 

Pappiz = 0.444) 

a 





-0.028+0.013 

-0.76+0.36 

b 



-0.032+0.038 

-0.85+1.02 

-0.025+0.013 

-0.67+0.36 

c 





0.118+0.013 

3.16+0.36 

d 

0.089±0.152 

2.38±4.07 

0.021+0.038 

0.57+1.02 

0.040+0.013 

1.06+0.36 

e 

0.194±0.152 

5.24+4.07 

0.409+0.057 

10.94+1.52 



f 

-0.072±0.152 

-1.93+4.07 

0.047+0.103 

1.25+2.75 



g 

-0.222±0.177 

-5.94+4.74 






Notes: (3 ap p(z = 0.444) is the corresponding apparent speed for /i at z=0.444. Considering the 
redshift range of 3C 66A (0.335 < z < 0.444), /3 app should be 0.18(3 app (z = 0.444) < /3 app < 
Pappiz = 0.444) 
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between 1.5 and 8 mas from the core. Following Lu et al. (2011), we could put constrains on basic 
jet parameters using the apparent speed, e.g. T min = ^1 + J3l pp and 6 max = 2 arctanfl l(i app ). For the 
fastest component, e (8.49 < f3 app < 10.94), we have, 8.55 < Y min < 10.99, 10.4° < 6 max < 13.4°, 
indicating that the Doppler boosting effect may play quite a large role in the very high core 
brightness temperature and detection of y-ray emission. 

Apparent inward motions have been detected in a number of individual objects (e.g. Mrk 
421, Niinuma et al. 2012, and OJ 287, Sawada-Satoh et al. 2013). But large sample studies have 
shown that such motions are relatively rare (Kellermann et al. 2004, Lister et al. 2013). Lister et 
al. (2013) reported that 17 out of 887 components (< 2%) show inward motions and that most 
of inward motions are slow (< 100 fias yr _1 ) and occur within ~ 1 mas from the core. We have 
detected inward motions in both the innermost (comp, a and b ) and outermost region (comp, g ) of 
the jet in 3C 66A. The motion of component g is probably due to the change of internal brightness 
distribution in this feature because it is relatively weak and extended. The reason for inward 
motions of the innermost components, as suggested by previous studies (e.g. Kellermann et al. 
2004), could be: (1) the motion of a newly emerged component which is still very close to the core 
and could not be separated by current resolution; (2) highly curved jet motions that cross the line 
of sight (LOS). The inward motions of the innermost components in 3C 66A were first reported 
by us at > 2cr confidence (Zhao et al. 2013) and were then confirmed by the MOJAVE results at 
15 GHz (Lister et al. 2013). The MOJAVE data cover mainly the time after our observation and 
we can cross identify most of our components with the 15 GHz ones (a, b, c, e, and f in our result 
with 13, 11, 9, 2, and 1 in MOJAVE results) by the extrapolation of the fitted motions. This means 
despite of a system shift of the component positions, which is caused by frequency-dependent 
core positions (i.e. the core-shift effect), the components at the inner part of the jet is continuously 
showing apparent-inward motions for ~ 10 years (2001 to 2011) and no new components were 
detected in this area even during and after y-ray and radio flares (e.g. in October 2008, Abdo et 
al. 2011). This result rules out both the two above explanations because otherwise for the former, 



- 19- 


the new-born component would show up after such a long time and for the later, the highly curved 
motions should turn over and cross the LOS again so the motions would change to outward in 
10 years. Here we propose another possible explanation: non-stationarity of the core due to 
opacity change. The non-stationarity of the core has been detected for some sources (e.g Mrk 421, 
Niinuma et al. in preparation, and 3C 454.3, Bartel et al. 2012) and the opacity change in 3C 66A 
is supported by the 2 epoch core-shift measurements in Pushkarev et al. (2012). Fig[4]shows the 
2D position plot of components a and b at 22 GHz. The average trend for the two components 
(bold arrow) are similar, which also indicate the change of the core position. The reason such a 
trend is not obvious in the other components is that these outer components were ejected much 
earlier than the inner components and according to the previous works the reported speeds were 
much faster (Jorstad et al. 2001, 2005). So their own motions dominate the overall motion. While 
for the innermost components, they are slower, even stable, so the core motion dominates the 
observed relative motions. Such a similarity could also be found in the MOJAVE results of this 
source. Note that the 3 third component in MOJAVE results is also likely inward moving, while 
it shows outward motions in 2001-2005 (Figure 3). This could be explained as this component 
decelerated around 2005 and began to show a similar motion as the inner 2 components. 

Future precise measurements of the core-shift with multi-frequency astrometric observations 
and detailed analysis of the variability of the core-shift and the possible correlation with flux, 
spectral variabilities would be important to test this explanation. 


5. Summary 

We studied the parsec-scale jet properties of the TEV blazar 3C 66A by using multi-epoch 
VLBI observations at 2.3, 8.4, and 22 GHz from 2004 to 2005. Our maps of 3C 66A show similar 
morphology as previous works, i.e. a core-jet structure with bendings. The northmost component, 
k is identified as the core based on its flat spectrum and high brightness. The other components 


DEC (mas) 
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Fig. 4.— 2D position plot of Component a (left) and b (right) at 22 GHz. The triangles, circles, 
and squares represent the results from C07, B05, and this work, respectively. The arrow shows the 
average trend of the component motion. (0,0) represents the position of the core. 
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are found to be weaker and show steeper spectra. 

We studied the proper motion of jet components in 3C 66A over a timerange of more than 4 
years and we found the kinematics of this source is very complicated. We detected superluminal 
motions for 3 components. We found the innermost components show inward motions. By 
combining with the results at 15 GHz, we have ruled out the possibility of new-born components 
or highly curved jet trajectories as reasons of the inward motion. We argue that the possible reason 
could be non-stationarity of the core. Further observations are needed to unveil the nature of this 
phenomenon. 
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